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ABSTRACT 

A point source observed 8 AU in projection from /^Pictoris in U (3.8 /am) imaging in 2003 has been 
recently presented as a planet candidate. Here we show the results of L'-band adaptive optics imaging 
obtained at Keck Observatory in 2008. We do not detect /3Picb beyond a limiting radius of 0.29'^ 
or 5.5 AU in projection, from the star. If /3Picb is an orbiting planet, then it has moved > 0.12^' 
(2.4 AU in projection) closer to the star in the five years separating the two epochs of observation. 
We examine the range of orbital parameters consistent with the observations, including likely bounds 
from the locations of previously inferred planetesimal belts. We find a family of low-eccentricity orbits 
with semimajor axes ~ 8-9 AU that are completely allowed, as well as a broad region of orbits with 
e < 0.2, a > 10 AU that are allowed if the apparent motion of the planet was towards the star in 2003. 
We compare this allowed space with predictions of the planetary orbital elements from the literature. 
Additionally, we show how similar observations in the next several years can further constrain the 
space of allowed orbits. Non-detections of the source through 2013 will exclude the interpretation of 
the candidate as a planet orbiting between the 6.4 and 16 AU planetesimal belts. 

Subject headings: circumstellar matter - planetary systems - stars: individual(HD 39060) 



1. INTRODUCTION 

A variety of observations support the hypothesis that 
A5V star P Pic hosts an orbiting planet. Smith fc Terrile 
(1984) noted that the dust distribution in the scattered- 
light disk should be centrally cleared or else significant 
extinction would be observed due to the edge-on disk ori- 
entation. A clearing within 100 AU radius was inferred 
by observations of scattered light and the rmal emission 
that probed closer to the star ( Golimowski et al.||1993j 
[Lagage & Pantin 1994; Kalas & Jewitt 1995). Through 
simulations. Rogues e t al. (1994) demonstrated that cen- 
tral clearings within debris disks are possible due to mean 
motion resonances of planets with Mp > 5 M^. Spec- 
troscopic observations revealed transient, redshifted ab- 
sorption features, which are interpreted as rapidly evap- 
orating comets initia lly located at a 4:1 resona nce with a 
hypothetical planet ( |Beust fc Morbideni|2000[ ). A signif- 
icant brightening and dimming of the pnotornetric light 
curve in 1981 is modeled as a planet transit event which 
could be consistent w ith an 8AU orbit if the p lanetary 
radius Rp > 0.2 (Lecavelier Des Etangs et al. 1997). 
Finally, an apparent warp m the inner disk midplane 
has been attributed to a planetary perturbation, imply- 
ing t hat the orbit is not coplanar with the extended dust 



disk ([Burrows et al.|19 95||Mouillet et al .|1997 
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The most recent evide nce for a pla net orbiting 0Fic 
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PA = 31.8° ±1.3°. They find the apparent bright- 
ness of the source, uil' — 11.2 ± 0.3 mag, is consis- 
tent with a 6-12 Mj planet at the stellar age (12^4 Myr; 
|Zuckerman et al.||2001^. Assuming a heliocentric dis- 
tance of 19.3 pc, the projected separation of /3Picb cor- 
responds to 7.93 ± 0.15 AU. Considering this as the min- 
imum semimajor axis of a planet circu larly orbiting an 
= 1.75 Mq star ( |Crifo et al.||1997[ ), then the period 
P > 16.9 yr. 

2. OBSERVATIONS & DATA REDUCTION 

2.1. Observations and Basic Calibration 

We observed /?Pic with the "narrow" (?^10 maspix"^) 
mode of the NIRC2 camera on the Keck II telescope on 
2008 December 2. Exposures consisted of 100 coadds 
of 53 ms integrations in correlated double sampling of a 
512x512 pixel subarray. We integrated 667.8s on-target 
in the U band (Aq = 3.8 /im, AA = 0.7 jam). The adap- 
tive optics (AO) loop was closed with P Pic serving as its 
own wave front reference. Despite the low elevation and 
partial vignetting by the lower edge of the dome shutter, 
we achieved moderate AO correction, with a diffraction- 
limited core and first Airy ring visible in individual ex- 
posures. The innermost portion of the PSF core satu- 
rated the detector. The average airmass (AM) was 3.1. 
Telluric and instrumental thermal backgrounds are sig- 
nificant in L' band. We alternated between pFic and 
dithered sky exposures (with equivalent readout config- 
uration) every 25 integrations for a total blank-sky inte- 
gration of 662.5 s. The field rotator was placed in "verti- 
cal angle" mode, which keeps the point-spread function 
(PSF) orientation fixed on the detector while the field 
rotates at the parallactic rate. The PFic exposure se- 
quence spanned 26.1° of parallactic rotation. 
We observed the star FS 13 {tul' = 10.10 ± 0.03; 



Leggett et al. 2003) to facilitate photometric calibra- 
tion. I'he low elevation of the /?Pic observations rel- 
ative to those of FS 13 (mean AM = 1.1) affects our 
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calibration strategy. We must account for the increase 
in atmospheric extinction as well as the degraded see- 
ing and AO correction of the /3Pic measurements. The 
vignetting during those exposures further reduced the 
throughput. Our goal is faint point source detection, so 
we focus on calibration in apertures having sizes simi- 
lar to the PSF core. We measure FWHM = 88.5 mas in 
the reduced FS 13 image, and adopt a nominal circular 
photometric aperture of this diameter for photometric 
calibration and subsequent point-source sensitivity anal- 
ysis. We derived the photometric zero point for this aper- 
ture, accounting for extinction, degraded Strehl ratio, 
and vignetting, by scaling the image of FS 13 to match 
/3Pic. We extracted a l^'xl^' box around each star and 
masked the saturated portion of /3 Pic's PSF core. We 
least-squares fit using the difference between the scaled 
image of FS 13 and P ic. Given m^' = 3.47 mag for 
pFic (Koornneef 1983), we derive a photometric zero 



in annuli around the star. Since speckles in AO-corrected 
images do not follow Gaussian statistics (e.g. Fitzger- 



point of 21.44 mag (defined at IDNs ^) for these aper- 
tures. This gives a scaling to monochromatic flux density 
of 6.610 X lO-'^JysDN-i ( [Tokunaga fc Vacca||2QQ5j) . 

We calibrated the detector orientation with 1.6 /im 
narrow-band observations of the ~ 6''-separation binary 
HD 56986. We found a 0.42° offset between the observed 
orientation and t he ephemeris compute d from the or- 
bital elements in |Hartkopf et al. (2001), and used this 
offset to correct the orientation of the /3 Pic exposures. 
Fo r scale, we a dopt 0.009963 arcsecpix"^ as measured 
by |Ghez et al. ( 2008| ; the value that we derived from the 
binary was 1% smaller. 



binary 
2.2. 



PSF Subtraction and Point-Source Detection 
Sensitivity 



We subtracted the stellar PSF usin g the angular dif - 
ferential imaging algorithm of Lafreniere et al | (|2007l). 
For each image in the exposure sequence, one constructs 
a stellar PSF reference using an ensemble of other im- 
ages in the sequence. The reference image is formed in 
patches, with the values in each patch determined by 
least-squares fitting a linear combination of images in 
the ensemble. In each patch, the selection of ensemble 
images is restricted to those in which the field-of-view has 
sufficiently rotated to prevent self-subtraction of point- 
source companions while maintaining fidelity in stellar- 
speckle suppression. After constructing a reference PSF 
for each image, the subtraction residuals are rotated to 
a common sky orientation and combined. 

We processed the images prior to PSF estimation. For 
stellar centroid estimation, we used a method similar to 
that outlined by Marois et al. (2006). After choosing one 
image as a reference and hign-pass filtering it, we ob- 
tained relative centroids by cross-correlating it with the 
other images. The reference image centroid was com- 
puted via cross-correlation with a version of itself rotated 
by 180°. After centroiding, we subtracted a radial profile 
from every image, using a robust mean in each annulus 
to avoid biases from diffraction spikes. Finally, we ap- 
plied a high-pass filter using a boxcar median, 30 pix on 
a side. 

Evaluation of PSF-subtraction performance requires 
one to both estimate the residual noise level and gauge 
the suppression of signals from point sources. We as- 
sessed residual noise by examining the distribution of 
fluxes from 88.5-mas-diameter circular apertures placed 



aid & Graham||2006 ), the distribution of residual tlux at 



a given location may not b e Gaussian even aft er com- 
bination of several images arois et al. 2008fc. With 
this caveat in mind, our confidence limits are based on 
the sample standard deviation of aperture photometry in 
each annulus, presuming Gaussian statistics. 

To assess the suppression of signal of faint sources in 
the data, we inserted artificial point sources in the images 
prior to processing. These sources' PSFs are determined 
by the normalized image of FS 13. We then applied the 
same boxcar filter as in the pre-processing step above, 
and scaled the sources to 11.2 mag (the candidate flux). 
In the outer regions of the image, sources are placed in a 
radial spoke pattern, with spokes separated by 60° and 
sources placed every 20 pixels. To prevent biases from 
crowding close to the star, sources at radius r are re- 
quired to maintain an arcwise separation rAO > 50 pix. 
We then measured the flux recovered by the algorithm at 
each radius with aperture photometry. Sources closer to 
the star undergo greater suppression. We least-squares 
fit an ad hoc function to the fraction of recovered flux as 
a function of radius. 



fir) 



a|l + 6[exp(r/c) - 1] ^| \ 



(1) 



This function provides a reasonable fit to the recov- 
ered artificial source fluxes from 0.2" -2.5". Here, a 
is the fraction of flux passed at large radii r, c is a 
transition radius where throughput drops, and b gov- 
erns the sharpness of the curve. To illustrate, at 

r = (0.3'^0.6'^1.0^^1.5'^2.0'^ 2.5^0. we find /(r) = 
(0.42,0.63,0.77,0.86,0.91,0.93), respectively. We then 
computed the point-source sensitivity as a function of 
radius after dividing the noise profile by the signal sup- 
pression function /(r). 

The subtraction algorithm has several parameters, 
which we tuned through trial and error to maximize the 
point- source sensitivity close to the star. ,Lafreniere et al.] 
(2007) discuss these parameters in detail We adopted 
hid = 88.5 mas for the resolution element width, A^^^ = 0.5 
for the minimum rotation for an exposure to be usable as 
a reference, and Na = 100 elements in the optimization 
area. The subtraction zone widths were dr = 1.5 x W 
for 0-125 pix, and b x W beyond. We adopted ^ = 1 for 
the optimization region aspect ratio. 



3. RESULTS & ANALYSIS 

Our final PSF-subtracted image is shown in Figure [l] 
The appare ntly point-like rriL' = 11.2 ± 0.3 mag object 
detected by Lagrange et al. (2009a) has not yet been con- 
firmed in the literature. In this section, we consider the 
implications of our U image for the presence of a bound 
companion. Before our analysis of companion orbits, we 
briefly consider the alternatives of a persistent feature in 
the circumstellar disk and an unbound background ob- 
ject. 



3.1. Expected Source Locations 

At the 200 3 source position (r = 411 ± 8 mas, P. A. = 
31.8° ± 1.3°; [Lagrange et al.' 2009a), we do not detect 
any point sources consistent with the flux measured in 
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Fig. 1. — Keck L^-band PSF-subtracted image of (3 Pic. North is up, east is left. The dotted circle indicates the extent of our 3-cr sensitivity 
limit {my = 11.5 mag; cf. Fig. [2|, and the region interior to this limit is shaded. No point sources are detected outside this region. A 
square northeast of the star marks the position of the planet candidate detected by Lagrange et al. (2009a) at r ~ 0.41^^, PA~ 32° in 2003 
Nov. A triangle to the lower left of the star indicates the expected position of the source in our 2UU8 Dec data if it were a background 
object, though this hypothesis cannot be ruled out since such a source would not be detectable in this region. The circle to the lower right 
is the expected position of the companion if it were on a circular orbit and located at maximum elongation in 2003 Nov, which again is 
not ruled out by the current data. Solid arrows indicate the orientation of the outer disk midplane (Kalas & Jewitt 1995), while dotted 
ones mark the position angle of the secondary disk inferred bv Golimowski et al., ( (2006| ). The inset shows output of the PSF-subtraction 
procedure after the insertion of artificial 11.2 mag sources ( ^2.2| . 

the discovery image. Figure [l] shows an elongated peak 
(r = 395 mas, P. A. = 38°) located near this discovery po- 
sition. We measure m^' = 12.3 mag in an 88.5 mas aper- 
ture, which is 3.1 X above the residual 1-a sample stan- 
dard deviation at that radius ( ^2.2[ ) and 1.1 mag fainter 
than the putative planet detected in 2003. We do not 
have sufficient confidence to assess this as a signal de- 
tection given the likely deviation from Gaussian statis- 
tics noted above. If the 2003 source was due to a local- 
ized density enhancement in the circumstellar disk, then 
it has not persisted to be detectable at the same peak 
brightness along the same line of sight in the 2008 data. 

We do not detect any point source at the position ex- 
pected if the Lagrange et al. (2009a) detection was due 
to an unlikely background object; however this position 
falls within our sensitivity limit and is therefore not sig- 
nificant. The sensitivity curve shown in Fig. [2] shows that 
we would likely detect an 11.2 mag object beyond 0.29'^ 
This limit is also shown as a dotted circle in Fig. [l] A 
background source would be detectable at this sensitivity 
beginning in 2010 Jan. 

Our image shows no 11.2-mag point sources outside the 
0.29^' limit. In Fig. [T] we mark the expected location if 
the companion was observed at maximum elongation in 
an edge-on circular orbit in 2003. This location is within 




Fig. 2. — Keck 3-cr L^-band sensitivity limit for ^^Pic. This curve 
was derived from the standard deviation of 88.5-mas-diameter 
aperture photometry, corrected for the radial dependence of the 
throughput of the roll-subtraction algorithm ( |2.2| . We sho w the 
detection of an 11.2 ± 0.3 mag source at r ~ U.4r^ by Lagr ange] 
et al. (2009a). For this range of source brightness, we show the 
corresponding projected separations, and adopt that of the dark- 
est vertical dotted line (r > 0.29^^) as an upper limit to the allowed 
projected separation of the source. 
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Fig. 3. — Map of orbital parameters allowed by the L' observa- 
tions assuming the candidate is a bound companion. At each posi- 
tion in (a, e), we show the fraction of orbital periastron longitudes 
/a; that are allowed given a detection in 2003 and our 2008 non- 
detection. The magenta curve shows the constraint rperi = 6.4 AU, 
while the green shows rap = 16 AU. The region of intersection be- 
low both of these curves is the range of orbits that are not allowed 
to cross the corresponding Okamoto et al. (2004) belt locations. 

our sensitivity limit, so we are unable to rule out a body 
in such an orbit. 

3.2. Orbital Constraints 

Our non-detection constrains other possible orbital 
configurations for the (3 Pic planet candi date. Given the 
detection at the best-estimate position of [Lagrange et al.| 
( |2QQ9a| ) in 2003, what orbital parameters are allowed 
with our exclusion of sources at r > 0.29'^ in the 2008 
epoch? For simplicity, we restrict our consideration to 
edge-on geometries {i = 90°). Atomic gas emission, re- 
solved in both space and velocity, shows quasi-Keplerian 
rotation with a redshifted NE ansa. We presume the 
candidate planet shares this orbital direction and define 
the longitude of the ascending node Q to correspond to 
the P.A. of the 2003 detection. The remaining free ele- 
mental parameters are semimajor axis a, eccentricity e, 
and argument of periastron co. If we define the direction 
of Q to be the x axis, we then have the y axis pointing 
away from the star along the line of sight and co being 
measured from x through y (i.e. argx + iy). We con- 
structed a grid in (a,e,a;) and computed the possible 
projected positions {x coordinates) of the source in the 
latter epoch. In general, there are two solutions for the 
anomaly of the 2003 epoch at each (a, e^uj) owing to the 
near and far points of intersection between an edge-on 
orbit and the line of sight. These solutions correspond 
to two possible positions in 2008. 

The range of orbital parameters allowed by the data 
is shown in Figure |3] The shading corresponds to the 
fraction of allowed orbits, f^^ over the full range of co. In 
the absence of prior information on co, this fraction cor- 
responds to probability of an allowed orbit given a and e. 
A key feature of the allowed parameter space is a family 
of low-eccentricity orbits with semimajor axes ~ 8-9 AU 
that are completely allowed by the observations (i.e. no 
constraints on a;). The completely disallowed region on 
the left edge of the graph (where f^j = 0) corresponds 



to the requirement that the apastron distance must be 
at least as great as the separation observed in 2003. On 
the other side of the Figure, there is a broad region of 
parameter space (e < 0.2, a > 10 AU) corresponding to 
wider orbits where the planet has a larger physical sep- 
aration but smaller apparent separation in the current 
data. Here, f^j = 50%, as half of the orbits (the near 
solutions) have the apparent separation increasing with 
time. 

We further constrain the allowable orbits by consid- 
ering the locations of planetesimal belts inferred from 
peaks in the crystalline silicate fraction obse rved with 
spatiahy re solved mid-infrared spectroscopy (Okamoto 



et 



at 

inferred 



20041). By requiring non-crossing orbits, a belt 
AU constrains the periastron distance, while the 
belt at 16 AU co nstrains apastron. As noted 
by Freistetter et al. (2007), these belt locations may be 
uncertain by a lew AU due to relatively coarse spatial 
sampling. We show these peri- and apastron constraints 
in Fig. [3j Together, they define a wedge-shaped region 
i n (g, e) space and requir e e < 0.4. 

'Boccaletti et al.| (2009) present non-detections > 0.5^' 
in the H band diUd^OA" in Kg in the 2004 Novem- 
ber. Due to the modeling uncertainties in translating 
these non-detections to L'-band equivalents, we do not 
incorporate these results into our orbital phase space con- 
straints. 

Previous works have made predictions for the orbits 
of putative planets, which we compare to our presently 
allowed space. Orbital predictions have been derived 
from the velocity distribution of transient gas absorption 
events hypothesized to arise from star-grazing comet-like 
bodies disrupted from their parent belt by mean-motion 
resonance interactions with a p erturbing planet (Beust 



et 



al. 
"tHe' 



Beust & Morbidelh 2000 ). For a planet 
these authors predict e > 0.6, 



1998 

m tiie 4:1 resonance, 
5AU ^ a ^ 20 AU, and a periastron longitude, mea- 
sured from the line of sight, between -80° and -60°. In 
the coordinate system defined here, this translates to 
— 170° ^ u; < —150°. Within these constraints, the cur- 
rent data do not place strong limitations on the range 
of allowable orbits. In this region, we find that orbits 
with apastra rap ^ 9.4 AU must correspond to the far 
solutions, i.e. those with a true anomaly between 0° and 
9 0° in the 2003 o bservation. 

Mou illet et al.| ([1997) have studied how a planetary 
perturber can mfiuence the apparent warp in the disk 
midplane. We note that the nature of this structure 
is unclear, as i t may be due to an inclined seco ndary 



disk d Colimowski et al. 
"caletti et al.| |20U9[) . As' 

mass range ol 6-1 3 



2006 however see also Boc- 



agrange et al. (2009a) note, a 
Mj estimated with atmosphere mod- 
els, couple d with the age- mass- semimajor- axis constraint 
of Mouillet et al.J_|r997^ , corresponds to a - 7.6-9.7 AU 
Our data in l^'ig. | 3] are largely consistent with this range 



Freistetter et al. p007 ) have numerically studied the pos- 
sible architectures olt he planetary system given the belt 
structure inferred by Okamoto et al. (2004) as well as 
the predictions for planetary perturbers that can account 
for falling evaporating bodies and the warp. The favored 
single-planet scenario has a 2 Mj planet at a = 12 AU 
with e < 0.1. The orbital parameters are consistent with 
the current data, though the mass of such a planet must 
be reconciled with the apparent brightness and age. 
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Fig. 4. — Map of orbital parameters for the candidate allowed 
by hypothetical future non-detections. As in Fig. |3] these panels 
depict the fuj at each (a, e); however in this case we presume non- 
detection at the same sensitivity each year through the date shown 
in each panel. The candidate detection of Lagrange et al. (2009a) 
can only be ruled out of this parameter range if non-detections 
continue through 2013. 



Most recently, Lagrange et al. ([2009bl) present non- 
detections of the candidate planet in L' and Kg data from 
2009 Jan and Feb, respectively. Restricting the compar- 
ison to L' ^ we achieve similar sensitivity levels. The data 
were obtained at similar epochs and there is a general 
concordance between the resulting conclusions, though 
we explicitly consider non-circular orbits . Bot h [Le cave- 



(2009) and Lagrange 



lier Pes Etan^s fc Vidal-Madjar 

et al." |2009b| ) consider whether tne candidate is consis- 
tent with causing the transit-like event observed in 1981. 
Again restricting consideration to the L' observations, 
both groups find consistency with orbits having semima- 
jor axes of ~ 8 or ~ 17 AU. As can be seen in Fig. [3) our 
non-detection is consistent with both families of orbits. 

Additional observations of /3Fic in the mid-term can 
further constrain the nature of the candidate. In Fig. [4] 
we examine the fraction of allowed orbits at each (a, e) in 
between the 6.4 and 16 AU planetesimal belts presuming 
no planet is detected outside of 0.29 in the next several 
years. For a non-detection in late 2009, the range of al- 
lowed orbits will be significantly constrained (top panel), 
with the elimination of longer period, near-circular orbits 
(with Tap ^ 13 AU) occurring the following year. The 
2003 detection can only be judged spurious if, at these 
sensitivity levels, non-detections continue through 2013. 



4. CONCLUSIONS 

We observed /3Fic in the L' band on 2008 Dec 2 and 
did not detect any point sources. In our data, the planet 
candidate observed in 2003 would have been detected if it 
were as close as 0.29^' from the star. Our sensitivity limit 
alone does not exclude the possibility that the source is 
an unlikely background object, nor can it rule out a com- 
panion on a circular orbit if it was observed at maximum 
elongation in 2003. By mapping the phase space of allow- 
able orbital parameters, and considering the likely loca- 
tion between two previously detected planetesimal belts, 
we determine that the planet candidate is spurious only 
if it is not detected in similar V observations conducted 
through 2013. 
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